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Abstract 


A summary  is  presented  of  a study  of  critical  fluctuations  using 
laser  light  scattering  f’-'hniques  ut.der  contract  N00014-67-A-0239-00I4 
and  N00014-76-C-0304  with  ONR.  Highlights  of  the  accomplishments  arc- 
experimental  verification  exponential  dveay  of  critical  fluctuations 
and  of  the  hypothesis  of  dynamical  scaling,  experimental  observation  -r 
an  effect  of  a viscosity  anomaly  on  the  decay  rate  of  the  fluctuations 
and  experimental  evidence  for  the  principle  of  universality  of  the  crit- 
ical correlation  function. 
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t.  Introduction 


f\  1 

This  report  summarizes  the  research  and  accomplishments  conducted  at  th'> 
University  of  Maryland  under  contract  N00014-6 7- A-02 39-0014  from  August  1970 
till  August  1975  and  under  contract  N00C14-76-C-0304  from  / ugust  1975  till 
August  1976  with  the  Office  of  Naval  Research. 

The  research  was  concerned  with  studies  of  fluctuations  in  aaterlals 
near  a critical  point  phase  transition  using  laser  light  scattering  techr.iuut 
The  research  was  conducted  under  the  supervision  of  Professors  C.O.  Alley  and 
J.V.  Senders  as  co-principal  investigators  with  Professor  R.F.  Chang  as  their 
aain  collaborator  on  this  project.  The  experiments  were  conducted  by  Professor 
Chat.g  in  the  quantum  electronics  laboratory  of  .‘rofessor  Alley.  Professor 
Sengers  supervised  the  interpretation  of  the  experimental  results.  During 
the  course  of  the  project  the  investigators  were  assisted  by  cwo  postdoctoral 
associates,  namely  Dr.  C.  jendjaballah  (1972-1973)  and  Dr.  A.J.  Bray  (1974-1976) 
and  by  three  graduate  students,  namely  P.H.  Keyes,  M.  Wigdor  and  H.  Burstyn. 

II.  Summary  of  Accomplishments 

Critical  point  phase  traur-itiocs  occur  in  a large  variety  of  dicfer«*  ‘ 
kinds  of  physical  systems  such  as  fervraagnets  near  the  Curie  point,  a-.,ui- 
ferronagnets  near  the  Neel  point,  alloys  exhibiting  order-disorder  transitions 
(e.g.  amnonium  chloride),  binary  liquid  mixtures  near  the  critical  point  of 
mixing  and  fluids  near  the  gas-liquid  critical  point.  These  diverse  physical 
systems  have  a number  of  phenomena  in  common  near  their  respective  critical 
points.  vrt»e  important  common  featui.  » is  the  existence  of  large  thermal 
: luctuutions  in  all  these  systems  near  the  critical  point. 
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The  main  purpcie  of  the  research  was  to  study  the  nature  of  these  critical 
fluctuations.  The  fluctuations  may  be  characterized  by  a correlation  length 
which  diverges  at  the  critical  point.  This  implies  that  in  the  vicinity  of 
a critical  point  the  correlation  length  becomes  much  larger  than  the  range 
of  intermolecular  forces.  Hence,  the  fluctuations  near  the  critical  point 
become  insensitive  to  the  details  of  thesfc  interactions.  It  is  believe-, 
therefore,  that  the  critical  fluctuation  behavior  in  the  various  different 
systems  h:’.s  an  essential  similarity,  referred  to  as  the  principle  of  uni- 
versality of  critical  behavior. 

Th<  existance  of  diverse  physical  systems  displaying  analogous  critical 
phenomena  provides  the  experimenter  with  the  opportunity  r,>  choose  both  the 
best  material  and  the  best  technique  for  studying  'pecific  aspects  of  crit- 
ical phenomena.  The  size  ot  the  fluctuations  in  magnetic  sotids  far  fr<  * 
the  critical  point  can  be  readily  measured  with  neutron  scattering  (such 
experiments  were  inducted  simultaneously  by  Professor  Minkiewicz  ani  are 
continued  by  Prof*.  .»r  Lynn  at  Zhe  University  of  Maryland).  Howeve  , this 
technique  lias  serious  limitations  very  close  t.o  the  critic-  *.  point  because 
the  fluctuations  extend  over  distances  much  larger  than  the  wavelengths 
rhe  neu'r m3.  Moreover,  very  close  to  the  critical  point  all  experiments  on 
solids  are  limited  by  lattice  strains  which  result  from  impurities,  vacancies, 
etc.  On  the  other  hand,  fluctuations  in  gases  and  liquids,  which  continuously 
"anneal"  themselves  car.  be  studied  accurately  by  laser  light  scattering 
techniques.  However,  experiments  neur  the  gas-liquid  critical  point  are 
severely  affected  by  multiple  scattering.  Hence,  experiments  in  binary 
liquids  near  the  critical  point,  where  multiple  scattering  can  be  minimized 
by  refractive  index  matching,  seemed  to  offer  the  best  opportunity  for 
studying  the  fluctuations  in  ue  close  vicinity  of  the  critical  point.  Hence, 
it  was  decided  to  measure  the  critical  fluctuations  in.  the  binary  liquid 
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3-methylpentane-nitroethane.  The  wisdom  of  this  choice  has  been  confirmed 
by  the  experimental  results  as  we  shall  see  below. 

Here  we  give  a brief  survey  of  the  results  of  the  i ^search.  For  the 
details  of  our  accomplishments  we  refer  f>  the  publications  and  reports 
listed  in  Section  111. 

When  the  experiments  were  initiated.  Professor  Korenman  made  a theot  i 
cal  analysis  of  what  type  of  information  could  be  obtained  fro*  the  applica- 
tion of  sophisticated  photon  counting  techniques  near  the  critical  point  of 
fluids.*  At  the  same  time  Professor  Sengers  made  an  analysis  of  the  nature 

of  the  transport  processes  and  the  decay  rate  of  the  fluctuations  near  the 

2 3 

critical  point  of  fluids  on  the  basis  of  existing  experimental  data.  ’ In 
particular,  he  found  that  the  decay  rate  data  could  indeed  be  described  in 
terms  of  the  modern  concepts  'or  treating  cril ical  point  phenomena  provided 
the  decay  rate  was  separated  in  an  anomalous  and  a regular  part. 

The  correlation  length,  determining  the  spatial  extent  of  the  critical 
fluctuations,  can  be  measured  by  determining  the  intensity  of  scattered  light 
as  a function  of  the  scattering  angle.  The  decay  rate  can  be  obtained  fror 
Che  frequency  distribution  of  the  scattered  light.  A theory  for  the  deca\ 
rate  of  the  critical  fluctuations  was  originally  developed  by  Kawasaki  at 
Tonplc  University  and  Ferrell  at  the  University  of  Maryland.  This  theory 
related  the  decay  rate  to  the  correlation  length  and  the  shear  viscosity  of 
the  system.  In  the  literature  it  was  claimed  that,  this  theory  gave  an  ade- 
quate account  of  the  experimentally  observed  decay  rate.  However,  due  to 
the  difficulty  of  obtaining  accurate  light  scattering  data,  and  hence  accurate 
data  for  the  correlating  length,  the  theory  was  always  applied  using  adjus- 
table parameters  for  the  temperature  dependence  of  the  correlation  length. 

In  our  opinion,  in  order  for  a crucial  test  of  the  theory,  it  was  imperative 
that  the  correlation  length  and  the  decay  rate  of  the  critical  fluctuations 
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both  be  measured  simultaneously  for  the  same  physical  system.  This  was  the 
purpose  of  our  first  series  of  experiments  completed  in  1972  and  reported  in 
references  4,5,6.  These  exp^iments  showed  for  the  first  time  that  the 
theoretical  expressions  derived  by  Kawasaki  and  Ferrell  in  fact  were  in- 
sufficient to  account  for  the  experimentally  observed  decay  rate  of  the 
critical  fluctuations.  This  result  was  caused,  in  part,  by  effects  from  a, 
anomalous  viscosity  behavior  on  the  decay  rate  and  it  led  to  considerable 
activity  by  a number  of  theoretical  investigators  at  the  University  of 
Maryland  as  well  as  elsewhere  and  several  refined  formulations  were  proposed 
yielding  improved  agreement  with  the  experimental  data.** 

These  theories  suggested  that  also  the  assumption  of  exponential  decay 
of  the  fluctuations  might  break  down  close  to  tK  :ritical  point.  With  an 
improved  photon  counf iug  detection  procedure  we  started  to  investigate  the 
nature  of  the  time  dependence  of  the  critical  fluctuations  more  accurately 
and  we  were  able  to  verify  the  exponential  nature  of  this  decay  up  to  very 
close  to  the  critical  point. ^ 

Nevertheless,  a number  of  questions  remained  unanswered.  First,  al- 
though the  new  tefined  theories  led  to  an  improved  agreement  between  theo. . 
and  experiment,  the  various  theories  (sometimes  referred  to  as  mode-coupling 
and  decoupled-mode  theories)  yielded,  in  fact,  predictions  different  from 
each  other.  However,  it  was  impossible  to  discriminate  between  these  theories 
based  on  the  existing  experimental  capability.  Secondly,  the  long  range 

character  of  the  critical  correlation  function  is  characterized  by  an  exponent 

If 

r which  appeared  to  be  a ve.  * elusive  property  to  be  measured  experimentally. 
Finally,  there  remained  the  nagging  question  whether  the  difference  between 
the  critical  exponents  deduced  experimentally  for  fluid  systems  and  the 


theoretical  values  calculated  for  the  Ising  model  indicated  a small  but 
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fundan*nta]  difference  between  fluid  systems  and  lattice  systems.  It  was 
eviden:  that  in  order  to  answer  these  questions  satisfactorily,  considerably 
more  accurate  light  scattering  intensity  measurements  would  be  required. 

This  g>al  has  been  the  object  of  our  research  during  the  past  three  years 
and  we  are  pleased  that  we  have  succeeded  in  obtaining  light  ^catt?rt;^; 
data  with  a precision  of  0.2Z  which  represents  an  improvement  of  the  pret 
by  one  order  of  magnitude. ^ 

Wien  data  are  to  be  obtained  with  such  precision,  many  small  effects, 
ptcvioisly  ignored,  must  be  taken  into  account.  Most  important  is  a precise 

,q 

assessaent  of  the  effect  of  double  scattering  which  was  studied  extensively . 

Tie  new  experimental  results  are  reported  in  the  Appendix  of  this  ro- 
't,  The  interpretation  of  these  results  was  gr.  tly  facilitated  by  the 
■ it  tint  Dr.  Bray  at  the  University  of  Maryland  had  been  able  to  develop  an 
’ urate  theoretical  formulation  for  the  correlation  function  of  the  Ising 
--del.  When  the  data  were  interpreted  in  terms  of  this  theoretical  corre- 
lation function,  it  appeared  that  the  critical  exponents  deduced  from  our  data 
for  the  binary  liquid  were  extremely  close  to  a recent  calculation  of 
" hese  exponents  for  the  tsing  model.  For  details  we  refer  to  the  Appeiwix, 

-.e  thus  have  obtained  at  the  conclusion  of  the  present  phase  of  the  project 
trong  experimental  evidence  that  the  critical  correlation  function  is  indeed 
iversal  and  the  same  in  fluids  and  solids.  This  is  the  first  tine  that 
eh  cxpeclmental  information  for  the  correlation  function  of  fluids  has  been 
4 talced.  We  may  confidently  expect  that  our  further  experimental  work  on 
: .t  ci itical  fluctiations  in  liquids  may  yield  valuable  information  for 
fluids  and  solid  systems  alike. 
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Abstract 


Using  light  scattering  techniques  we  have  measured  the  correlation 
function  of  a binary  liquid  near  the  critical  point  with  a precision  of 
0.2Z.  The  data  cover  a range  of  0.18  < k£  < 26  in  terms  of  the  scaling 
variable  k£  and  yield  the  critical  exponent  values  y * 1-240  ± 0.007, 
v * 0.625  ± 0.003  and  n = 0.016  ± 0.007.  The  results  indicate  that  binary 
liquids  belong  to  the  same  universality  class  as  the  Ising  model. 
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Hear  a second  order  phase  transition  the  intensity  of  scattered  radiation 
for  scattering  vector  1c  is  proportional  to  the  Fourier  transform  xOO  of  the 

A. 

order  parameter  correlation  function.  Sufficiently  close  to  the  critical 
reapers ture  T£  one  expects  the  scaling  fora 

*(k)  * r t"Tg(lc>)  (1) 

Where  t * (T-T  )/T  , y the  susceptibility  exponent  , 5 the  correlation  length 
c c 

and  g(x)  the  correlation  scaling  function^. 

According  to  the  universality  hypothesis,  systems  having  the  saae  basic 

syaaetries  are  expected  to  have  identical  critical  exponents  and  scaling 

2 

functions,  and  are  said  to  belong  to  the  sane  universality  class.  It  is 
widely  assumed  that  universality  classes  for  homogeneous,  isotropic  systess 
may  be  assigned  according  to  the  spatial  dimensionality  d and  the  number  of 
components  n of  the  order  parameter.  This  assumption  implies  that  fluids  near 

the  gas-liquid  critical  point  and  binary  liquids  near  the  consolute  point  ; 

l 

should  belong  to  the  same  universality  class  (d  * 3,  n * 1)  as  the  three-  ; 

I 

dimensional  Ising  model.  The  critical  exponents  for  this  universality  clas* 

3 k 

have  been  calculate  1 using  high-teaperature  series  expansion  techniques  * 
and  renormalization  group  techniques'* *^.  The  various  calculations  yield 
numerically  similar  results  although  some  small  unresolved  discrepancies 
do  exist. 

During  the  past  years  the  values  reported  for  the  critical  exponents  of 
7 3 

gases  as  well  as  binary  liquids  appeared  to  deviate  from  those  attributed  to 
the  Ising  model.  However,  recent  experiments  of  Hoeken  and  Moldover  on  Xe, 

SF.  and  CO.  extremely  close  to  the  gas-iiquid  critical  point  indicated  that 

t>  L 
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the  exponents  do  approach  Ising  like  values  . In  addition,  recent  co- 
existence curve  data  obtained  by  Greer  for  the  isobutyric  acid  and  water 
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mixture  showed  good  agreement  between  the  experimental  and  the  Ising  value 

r 10 

for  the  exponent  g . 

In  order  to  further  investigate  the  hypothesis  of  universality,  with 
particular  emphasis  on  that  for  the  critical  correlation  function,  we  have 
obtained  a series  of  accurate  light  scattering  measurements  for  the  binary 
liquid  3-methylpentane-nitroethane.  The  measurements  enable  us  to  determine 
the  critical  exponents  and  the  correlation  scaling  function-  The  values 
obtained  for  the  critical  exponents  are  very  close  to  those  obtained  from 
the  Ginzburg-iandau-Wilson  formulation  of  the  Ising  model  by  Baker  et  al.^ 
indicating  that  binary  liquid  mixtures  do  indeed  belong  to  the  same  univer- 
sality class  as  the  Ising  model. 

Using  photon  counting  techniques  the  scattered  intensity  at  a scattering 

angle  of  90°  was  measured  as  a function  of  temperature.  Data  points  taken  in 

—6  .3 

the  range  1 x 10  < t < 3 x 10  , corresponding  to  0.18  < k£  < 26,  were  used 

in  the  analysis.  In  order  to  obtain  accurate  absolute  scattering  data  each 
scattering  intensity  reading  was  normalized  with  respect  to  the  incident 
intensity  measured  with  the  same  photodetector.  The  r.m.s.  uncertainty  of 


each  data  point  is  roughly  0.2%,  determined  from  twenty  consecutive  readings 


collected  at  100-second  counting  intervals  at  a fixed  temperature.  The  tem- 
perature was  stabilized  to  vir'nin  0.2  millidegrees.  The  concentration  of 


the  mixture  was  within  0.4%  of  its  critical  value.  Experimental  details  will 
be  presented  elsewhere. 

When  the  precision  of  the  data  is  a fraction  of  one  percent,  many  smell 
effects,  previously  ignored,  must  be  taken  into  account.  Therefore,  the 
data  are  corrected  for  the  following  effects:  (1)  turbidity,  (2)  double 

scattering**,  (3)  temperature  dependence  of  the  intensity  prefactor  TP  *(3c/3c)‘ 


where  p is  the  density,  c the  dielectric  constant  and  c the  concentration 
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(order  paraaeter) , and  (4)  entropy  fluctuations,  Brillouin  scattering  and 

cell  wall  scattering.  These  effects  were  either  aeasured  or  estimated  from 

known  physical  properties  of  the  mixture.  The  magnitudes  of  the  turbidity  and 

double  scattering  corrections  are,  respectively,  about  42  and  12  at  « 10, 

22  and  0.42  at  k£  * 3,  and  become  progressively  smaller  as  decreases.  The 

tcaperature  dependence  of  the  prefactor  causes  it  to  deviate  from  its  value  at 

T^  by  0.12  at  k£  * 1 and  22  at  k£  * 0.18.  The  extraneous  contributions  from 

entropy  fluctuations  and  Brillouin  scattering  are  estimated  to  be  roughly  two- 

thirds  of  the  scattering  intensity  at  T-T^  * 10  K and  are  nearly  constant  in 

the  temperature  range  of  the  experiment. 

Past  determinations  of  the  exponent  n from  scatte-ing  experiments  have 

been  difficult,  partly  due  to  the  lack  of  a scaling  function  g(x)  of  suffic- 
13 

ient  accuracy  . This  deficiency  may  be  remedied  by  using  a "truncated  Fisher- 

14 

Langer”  approximant  recently  proposed  by  one  of  the  authors  . This  repro- 
duces to  high  accuracy  the  theoretically  known  Ising  model  correlation  functions 
in  2 and  (4-e)  dimensions.  To  facilitate  the  analysis  we  use  a form  linearized 
in  n.  A numerical  check  confirms  that  the  linearization  is  valid  for  the  smaix 
values  of  n considered  here.  The  reciprocal  of  the  scaling  function  is  given  by: 


-n/2 


g'^x)  - 1 + x2(l  + f) 


where  s^Cx)  is  given  by 


-V*> ' t 


du 


J 


u(u2  + x2/9) 


[1  - F(u)] 


(2) 


(3) 
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Here  F(u)  is  the  spectral  function  ' * , truncated  so  as  to  vanish  for  u < 2, 
which  is  derived  from  the  Fisher-Langer  scaling  function*^ 


4 


^(x)  - (Cj^/x2”15)  (1  + C2/x(1  ' a)/v  + C3/x1/V),  for  x »1, 

through  the  relation  la  g^(ijxj)  * [sin^n^J/C^]  jx!2*^  F({x|/3)..  The 
spectral  function  may  b?  "fins— tuned"  by  adding  a constant  contribution,  of 
strength  w i 0 in  the  range  l<u<2.  This  contribution  leads  to  the  final  term 
in  the  square  brackets  in  Eq.  (2), 


The  data  analysis  is  carried  out  by  means  of  a non-linear  least  squares 
fitting  procedure  using  five  free  paraaeters:  an  overall  proportionality 

constant  Iq.  the  critical  exponents  v and  n,  the  correlation  length  amplitude 
and  an  extraneous  intensity  contribution  tl.  The  best  fit  -values  for 
AI  turn  out  to  be  very  close  to  the  estimates  discus sc ! above.  Preliminary 
analysis  indicated  that  the  best  fit  values  of  n and  v are  cl t-oe  to  0.015  and 
0.625  respectively.  Subsequently  for  the  final  analysis  we  used  an  s^(x) 
evaluated  with  v * 5/8  and  a * 1/8,  the  latter  derived  from  t;,e  hyperscaling  re- 
lation 2 - a « 3v,  and  the  values  ?j  * 1/54  = 0.0185,  C2  - 1.773  and  C3  « -2.745 

2 17 

derived  fron  the  e-expansion  to  0(e  ) . The  use  of  these  values  leads  to 

14 

a spectral  function  of  the  expected  shape  . Since  the  best  fit  values  of 
the  critical  exponents  were  found  to  be  insensitive  to  the  choice  of  v within 
a physically  reasonable  range  (0  $ v < F(2)  ^ 0.17),  the  value  w*»0  was  used 
in  the  final  analysis.  The  closeness  of  the  best  fit  values  of  v and  n to 
those  used  in  calculating  s2<x)  supports  a posteriori  the  self-consistency 
of  the  analysis. 


Eight  experimental  runs  were  available  for  this  analysis,  the  results 
of  which  are  presented  in  Table  I.  The  values  for  y and  n. , which  are  the 
basic  exponents  calculated  by  Baker  et_  al.^  were  deduced  from  v and  n using 
the  scaling  law  y - v(2  - n)  ar  * the  definition  n * q - 2 + 1/v.  The  first 
six  runs  were  obtained  with  the  scattering  beam  located  at  the  level  where  the 
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meniscus  would  appear.  In  the  last  two  runs  this  position  was  varied  by  2 and 
5 m'  from  the  earlier  position.  The  results  appear  to  be  independent  of  the 
level  of  the  beam  confirming  the  absence  of  gravity  effects.  The  average 
values  and  r.a.s.  errors  of  the  critical  exponents  from  the  eight  runs  are 

v - 0.625  ± 0.003  y « 1.240  ± 0.007  (4) 

n « 0.016  ± 0.007  tu,  —0.384  ± 0.010 

A 6 18 

When  an  extended  scaling  correction  factor  (1  + Ct  *)  with  Aj  * 0.5  * 

was  included  in  the  analysis,  we  obtained  C » 0.3  - 0.3  with  insignificant 
shifts  in  the  exponent  values  indicating  that  the  data  are  within  the  range 
of  asymptotic  simple  scaling. 

Our  results  are  rather  similar  to  the  exponent  values  calculated  by 
baker  et  al.  from  the  Ginzburg-Landau-Wilson  model  using  the  Callan-Symanzik 
equation*’ : 

v « 0.627  ± 0.01  y - 1.2410  t 0.002  (5) 

n - 0.021  ± 0.02  Hu  “ -0.3843  ± 0.003 

The  data  were  also  compared  with  the  predictions  of  high  temperature  series 

# 

expansions  * by  fixing  y and  n at  0.638  and  0.041  and  using  an  appropriate 
14 

scaling  function  . For  c caparison,  two  deviation  plots  of  a representative 
experimental  r in  (run  no.  3)  are  presented  in  Fig.  1.  Whereas  the  upper  plot 
with  v * 0.625  and  n * 0.01"  does  not  show  any  significant  deviations,  small  but 
noticeable  systematic  deviations  do  appear  in  the  lower  plot  with  v * 0.638  and 
rj  * 0.041.  Our  dat£  thus  seems  to  support  the  Ising  exponents  as  computed  by 
Baker  et  . 

Our  exponent  y is  also  in  good  agreement  with  the  value  of  this  exponent 

found  by  Hocxen  and  Holdover  for  gases  asymptotically  near  the  critical 
o 

point  . However,  unlike  gases  where  the  asymptotic  behavior  is  not  approached 
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^ unless  t * 10  we  find  that  for  binary  liquids  the  range  of  asywptotic 

1 .3 

m behavior  extends  to  t * 10  in  agreement  with  Greer's  analysis  of  the  co- 

I existence  curve  of  a binary  liquid***. 

In  conclusion  we  note  that  our  analysts  daeimines  both  the  exponents 

I and  the  correlation  scaling  function.  The  close  agreement  between  our  results  - 

■ j 

_ and  those  of  Baker  et  al.  indicates  that  binary  liquids  and  the  Ising  Model 

» do  belong  to  the  sane  universality  class.  Moreover,  unliae  gases  this 

■ universal  behavior  is  observable  in  a temperature  range  that  is  readily 
accessible  experimentally . 
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TABLE  I.  EXPERIMENTAL  CORRELATION  FUNCTION  PARAMETERS 
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Fig.  1 Percentage  deviations  bei»«en  the  measured  scattering  intensity 
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